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TABLE. 1 . Model parameters. 


Models 

r, mm 

Y.s-> 

D 

T c , K 

T bl .K 

1 

1 

10-15 

D AI,Op* 

1050 

1321 

2 

10 

10-15 

D Al,Opx 

1050 

1321 

3 

1 

10-15 

^Fc.Gt 

1050 

1321 

4 

1 

10-16 

D Al,Op* 

1050 

1321 

5 

10 

10-16 

D AJ,Opx 

1050 

1321 

6 

10 

10-16 

I^Al,Opx 

1350 

1521 

7 

1 

10-16 

D Al ( Opx 

1150 

1521 

8 

10 

10-16 

D Fc,Gt 

1050 

1321 


crustal -heat-production dominated. For y = 1(H 5 s -1 , the observed 
maximum elevation of mountain belts can be explained as a 
consequence of disequilibrium phase boundary depth for a wide 
range of physical parameters, although a comparatively young age 
for Maxwell Montes (50 m.y.) is implied. For the lesser horizontal 
strain rate of 10~ 16 s -1 , only limited parameter values for thermal 
models are allowed. 
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RESULTS OF A ZONALLY TRUNCATED THREE-DIMEN- 
SIONALMODEL OF THE VENUS MIDDLE ATMOSPHERE. 
M. Newman, CIRES, University of Colorado, Boulder CO 80309, 
USA. 

Although the equatorial rotational speed of the solid surface of 
Venus is only 4 m s -1 , the atmospheric rotational speed reaches a 
maximum of approximately 100 m $~ l near the equatorial cloud top 
level (65 to 70 km). This phenomenon, known as superrotation, is 
the central dynamical problem of the Venus atmosphere. We report 
here the results of numerical simulations aimed at clarifying the 
mechanism for maintaining the equatorial cloud top rotation. 

Maintenance of an equatorial rotational speed maximum above 
the surface requires waves or eddies that systematically transport 
angular momentum against its zonal mean gradient. The zonally 
symmetric Hadley circulation is driven thermally and acts to reduce 
the rotational speed at the equatorial cloud top level; thus wave or 
eddy transport must counter this tendency as well as friction. 
Planetary waves arising from horizontal shear instability of the 
zonal flow (barotropic instability) could maintain the equatorial 
rotation by transporting angular momentum horizontally from 
midlatitudes toward the equator. Alternatively, vertically propagat- 
ing waves could provide the required momentum source. The 


relative motion between the rotating atmosphere and the pattern of 
solar heating, which has a maximum where solar radiation is 
absorbed near the cloud tops, drives diurnal and semidiurnal ther- 
mal tides that propagate vertically away from the cloud top level. 
The effect of this wave propagation is to transport momentum 
toward the cloud top level at low latitudes and accelerate the mean 
zonal flow there. 

We employ a semispectral primitive equation model with a zonal 
mean flow and zonal wavenumbers 1 and 2 [1]. These waves 
correspond to the diurnal and semidiurnal tides, but they can also be 
excited by barotropic or baroclinic instability. Waves of higher 
wavenumbers and interactions between the waves are neglected. 
Symmetry about the equator is assumed, so the model applies to one 
hemisphere and covers the altitude range 30 to 1 10 km. Horizontal 
resolution is 1.5° latitude, and vertical resolution is 1.5 km. Solar 
and thermal infrared heating, based on Venus observations and 
calculations drive the model flow [2]. Dissipation is accomplished 
mainly by Rayleigh friction, chosen to produce strong dissipation 
above 85 km in order to absorb upward propagating waves and limit 
extreme flow velocities there, yet to give very weak Rayleigh 
friction below 70 km; results in the cloud layer do not appear to be 
sensitive to the Rayleigh friction. The model also has weak vertical 
diffusion, and very weak horizontal diffusion, which has a smooth- 
ing effect on the flow only at the two grid points nearest the pole. 

Simulations were carried out with uniform background angular 
velocity equivalent to an equatorial speed of Uq, where Uq was varied 
between 50 and 75 m s -1 . Row with this angular velocity was the 
initial condition for half of the simulations. The initial condition for 
the other half was obtained by adding to this background rotation a 
horizontally uniform, cyclostrophic-balanced component with zero 
additional zonal velocity at 30 and 1 10 km and a smooth increase to 
a maximum addition of 50 ms - * at 65 km on the equator. Model runs 
were also carried out in which the coefficient of vertical diffusion 
u was varied. Cases were run for 350 simulated (Earth) days, by 
which lime a statistically steady state was reached. We present 
averages for the last 40 days of each run. 

In the resulting mean zonal flow, the equilibrated equatorial 
wind maximum was typically between 90 and 105 m s _1 , and a jet 
developed near 40° latitude. The tides, particularly the semidiurnal 
tide, acted to balance (upgradient/downgradient) vertical advection 
by the Hadley cell updraft (below/abovc) the low -latitude zonal 
wind maximum. Experiments in which Uq was vaned indicated that 
the shape of this vertical jet (i.e., the vertical wind shear) is less 
sensitive to the background rotation than is the value of the speed 
maximum. This suggests that any theory that describes the role of 
the thermal tides in maintaining the equatorial rotational wind 
structure against vertical advection cannot oversimplify the vertical 
wind structure. Vertical diffusion acted to counter the tidal accelera- 
tion at cloud top, producing a slower zonal wind speed. For example, 
for two runs employing the sheared 60 m s _l background rotation, 
reduction of u from 2.5 to 1.0 m 2 s~* produced an increase of more 
than 10ms _I in the equatorial zonal wind maximum. The diurnal 
tide transported angular momentum horizontally from the region of 
the midlatitude jet toward lower latitudes, acting to smooth the 
zonal wind profile between the mid latitude jets and the equator; 
thus, the jet does not become significantly barotropic ally unstable. 

The tides also acted to weaken the Hadley circulation through 
both their meridional and upward heat fluxes. At the equator, the 
vertical convergence of the upward heat flux compensates part of 
the zonally averaged solar heating in the cloud top region. This 
effect reduced the mean equatorial cloud -level updraft by half, as 
compared to zonally symmetric model runs. Thus the tides act to 
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Fig. 1. Comparison of model distributions of cloud top zonal and meridional 
winds with local time, (a) and (b), with distributions derived from cloud drift 
measurements, (c) and (d) [3}. The heavy black lines mark the region of UV 
cloud observations. Model results are averages from the last 40 days of the 
350-day simulation with initial sheared angular velocity, for u 0 = 60 m s“ l and 
v = 2.5 m 2 s* 1 . 

maintain the mean zonal flow at cloud top level in two ways: directly 
through tidal fluxes of angular momentum, and indirectly, by 
reducing the thermal forcing of the Hadley cell. 

At the cloud top level, the model zonal wind has a distinct 
minimum during daytime such that the midlatirude jet is substan- 
tially diminished from the zonal mean value. The daytime meridional 
wind, on the other hand, is strongly poleward and is at a maximum 
in early afternoon (Fig. 1). This is a robust result that is consistent 
with both the near absence of midlatitude jets and the strong 
meridional winds deduced from cloud drifts. Cloud drifts can be 
measured only during the period from about 0800 to about 1600 
hours local time, but they have been used to infer zonally averaged 
zonal and meridional winds [3]. As shown by the results, this is 
inappropriate and would lead to spurious estimates. 

Also present are nonlidal wave components: a wavenumber 2 
component with typically a 2.6-day period, and a weaker wavenumber 
1 component with the same angular phase speed and 5.2-day period. 
These are confined mainly to midlatitudes and are due to baroclinic 
and barotropic instabilities of the zonal mean flow in the lower cloud 
region. These features may be related to a wavenumber 1 5.2-day 
spectral feature that has been detected in cloud-level observations. 
The wavenumber 2 feature became stronger relative to wavenumber 
1 as the background rotation was increased between runs, which 
caused the jet to be centered at slightly lower latitude. This relative 
wavenumber 2 dominance may also be a consequence of the 
omission of wave-wave interactions. These waves did not contrib- 
ute to maintenance of the zonal mean flow near the equator, but in 
mid to high latitudes they did act to transport momentum from the 
middle to lower cloud region, and thus possibly contributed to the 
slight (-5 K) temperature rise in the cloud-top polar region. Al- 
though the model did not duplicate the observed cloud top polar 
warm region, it appears that a higher-latitude jet in themiddle cloud 
region could act to both increase the wavenumber 1 component and 
duct the waves more toward the pole, inducing warmer temperatures 
there. 

To determine the robustness of these results, further model runs 
are currently being conducted to more completely explore the 
parameter space described above. These results also suggest that a 
Venus GCM will either need sufficient resolution to capture the 


tides, or else a parameterization scheme that can adequately capture 
their effects. 
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A MODEL FOR THE FORMATION OF THE EXTENDED 
PARABOLOIDAL HALOS AROUND SOME IMPACT CRA- 
TERS ON VENUS. W. I. Newman 1 , E. M. Jones 2 , D. B. Campbell*, 
and N. J. S. Stacy 3 , Apartments of Earth and Space Sciences, 
Astronomy, and Mathematics, University of California, Los Ange- 
les CA 90024-1567, USA, 2 Earth and Environmental Sciences 
Division, Los Alamos National Laboratory, Los Alamos NM 87545, 
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Many parabolic- shaped extended impact crater-related features 
have been found in Magellan synthetic aperture radar and emissivity 
data covering much of the surface of Venus. They are oriented east- 
west with the apex to the east and the impact crater located just west 
of the apex. A model for the formation of the parabolic features is 
developed based on the injection of small particles into the upper 
atmosphere at the time of impact, and their subsequent transport to 
the west by the east- west zonal winds. Fallout times from 50 km in 
the Venus atmosphere for particles of this size are about two hours, 
allowing westerly drifts of several hundred kilometers for zonal 
winds of 50 to 100 m s* 1 ^ - - 
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VENUS: PRELIMINARY GEOLOGIC MAPPING OF 
NORTHERN A TLA REGIO. A. M. Nikishin 1 andG. A. Burba 2 , 
l Geological Faculty, Moscow State University, 119899, Moscow, 
Russia, 2 Vemadsky Institute, Russian Academy of Science, 1 1 7975 , 
Moscow, Russia. 

A preliminary geologic map of C 1 sheet 1 5N 1 97 was compiled 
according to Magellan data (Fig. 1). Northern Atla Regio is domi- 
nantly a volcanic plain with numerous volcanic features: radar- 
bright and -dark flows and spots, shield volcanos, volcanic domes 
and hills with vaned morphology, and coronalike constructions. 
Tesserae are the oldest terrains semiflooded by plain materials. 
There are many lineated terrains on this territory. They are inter- 
preted as old, partly buried ridge belts. Lineated terrains have 
intermediate age between young plains and old tesserae. Ozza Mons 
and Sapas Mons are the high shield volcanos. 

The prominent structure of northern Atla Regio is Ganis Chasma 
rift. The rift dissected the volcanic plain and evolved nearly contem- 
poraneously with Ozza Mons shield volcano. Ganis Chasma rift 
valley is highly fractured and bounded by fault scarps. There are a 
few relatively young volcanic features in the rift valley. The rift 
originated due to 5-10% crustal extension and crustal subsidence 
according to analysis of fracturing and rift valley geometry, Ganis 
Chasma is characterized by rift shoulder uplifts. Geological struc- 
tures of Atla Regio and Beta Regio are very similar as assumed 
earlier [1,2]. 
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